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Summary 

Pseudomonas putida DOT-T1E is an organic solvent 
tolerant strain capable of degrading aromatic hydro- 
carbons. Here we report the D0T-T1E genomic 
sequence (6 394 153 bp) and its metabolic atlas 
based on the classification of enzyme activities. The 
genome encodes for at least 1751 enzymatic reac- 
tions that account for the known pattern of C, N, P 
and S utilization by this strain. Based on the poten- 
tial of this strain to thrive in the presence of organic 
solvents and the subclasses of enzymes encoded in 
the genome, its metabolic map can be drawn and a 
number of potential biotransformation reactions can 
be deduced. This information may prove useful for 
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adapting desired reactions to create value-added 
products. This bioengineering potential may be real- 
ized via direct transformation of substrates, or may 
require genetic engineering to block an existing 
pathway, or to re-organize operons and genes, as 
well as possibly requiring the recruitment of 
enzymes from other sources to achieve the desired 
transformation. 

Bacteria of the genus Pseudomonas are motile Gram- 
negative bacteria characterized by high metabolic versa- 
tility, and aerobic respiration, although a few strains of 
different species are able to use nitrate as a final electron 
acceptor (Palleroni, 2010). Pseudomonads are ubiquitous 
soil and water microorganisms that colonize many differ- 
ent environments and, consequently, have diverse life- 
styles. Strains of the species Pseudomonas putida are 
frequently soil inhabitants and are important in organic 
matter recycling in nature; they have a high bioremedia- 
tion potential because they often carry genes to deal with 
natural and xenobiotic chemicals (Nelson era/., 2002; 
Caballero era/., 2005; van Dillewijn etal, 2007; Arias 
etal., 2008; Segura etal., 2009a,b). The key to the ubiq- 
uitous distribution of these bacteria is not only their meta- 
bolic potential, but also the arsenal of regulatory genes 
that allow them to adapt to changes in the environment 
(Sashidhar and Podile, 2009; Wu etal., 2011). A few 
P. putida strains, namely, S12, Idaho and DOT-T1E 
(Weber etal., 1993; Ramos etal., 1995; Pinkart etal., 
1996), are able to thrive in the presence of toxic solvents 
(e.g. decanol, octanol, toluene, styrene), and these 
strains are considered extremophile microorganisms with 
great potential in bioremediation and in biocatalysis in 
biphasic systems (Ramos etal., 1995; Isken and de Bont, 
1996; Molina etal., 2011; Ramos etal., 2011; Tao etal., 
2011; Udaondo etal., 2012). Organic solvents are toxic to 
most microorganisms because they dissolve in the cell 
membranes, disorganize their structures and impair vital 
functions such as respiration, and the collapse in energy 
generation lead to cell death (Sikkema etal., 1995; 
Ramos etal, 2002; 2011). Solvent tolerance in P. 
putida DOT-T1E is a multifactorial trait that involves 
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chromosomal and plasmid encoded functions (Ramos 
etal, 2002; Segura etal., 2005; 2009b; Rodn'guez-Herva 
etal., 2007; Garcia etal., 2010; Molina etal., 2011). The 
first barrier to solvents involves a reduction in the perme- 
ability of the cell membrane via a fast cis to trans isomeri- 
zation of unsaturated fatty acids followed by a slower 
modification of phospholipid head groups (Keweloh and 
Heipieper, 1996; Junker and Ramos, 1999; Heipieper 
etal., 2001; 2003; Bernal etal., 2007; Pini etal., 2009). 
However, this reduction in permeability does not prevent 
entry of the solvents, which results in unfolding of proteins 
and the consequential function of a number of chaper- 
ones (Segura etal., 2005; Domi'nguez-Cuevas etal., 
2006; Volkers etal., 2006). The main mechanism under- 
lying solvent tolerance lies in the action of RND 
(resistance-nodulation-cell division) efflux pumps 
encoded on the host chromosome and on the pGRT1 
plasmid (Kieboom etal., 1998; Kim etal., 1998; Ramos 
etal., 1998; Mosqueda and Ramos, 2000; Rojas etal., 
2001; Rodn'guez-Herva etal., 2007; Godoy etal., 2010; 
Udaondo etal., 2012). Resistance to solvents is also 
modulated by the action of chromosomally encoded ABC 
efflux transporters that use energy to remove solvents 
from the cells to the outer medium (Kim etal., 1998; 
Garcia etal., 2010). 

Here, we present, the genome of the solvent tolerant 
P. putida DOT-T1 E strain obtained using the 454 technol- 
ogy. This microorganism uses a wide range of carbon, 
nitrogen, sulfur and phosphorous sources due to its wide 
metabolic potential. In addition we summarize previous 
knowledge on the biotransformation potential of this strain 
and how the properties and genomic information can be 
used to design new biotechnological processes. 

Genome sequencing, assembly, annotation and 
bioinformatic analysis 

The complete genome sequence of P. putida DOT-T1 E 
(GenBank accesion number CP003734) was determined 
by applying a strategy combining whole-genome-shotgun 
454-pyrosequencing on the genome sequencer FLX plat- 
form (20x coverage with 305 contigs, 257 354 nt biggest 
contig) and Sanger sequencing of PCR amplicons cover- 
ing gaps between contigs. In addition, after a first round of 
annotation, regions of lower quality as well as regions with 
putative frameshifts were resequenced from PCR ampli- 
fication of the dubious regions and the complete genome 
sequence was established. The genome of P. putida 
DOT-T1E has two circular replicons; a single chromo- 
some of 6 260 702 base pairs with a G+C content of 
61% (Fig. 1) and another of 133 451 base pairs corre- 
sponding to a plasmid previously named pGRT1 (Molina 
etal., 2011; Genebank HM626202.1, NCBI Reference 
Sequence: NC_015855.1). 



Using a combination of the Glimmer 3.03 software 
(Salzberg etal., 1998; Delcher etal., 1999), blast analy- 
sis and manual curation a total of 5803 ORFs were pre- 
dicted and annotated in the chromosome, of which 170 
have no significant homology (at E-value of < 10~ 5 ) to 
any ORF present among the sequenced Pseudomonas 
genomes. The analysis of ORFs of pGRT1 was performed 
previously and revealed it encodes for 126 proteins 
(Molina etal., 2011). In total the genome of DOT-T1E 
encodes 5721 proteins and 82 RNAs of which 58 corre- 
sponded to tRNAs. We analysed the GC skews of the T1 E 
chromosome, which is defined as the value of [G-C]/ 
[G+C] where G and C represent the local base frequen- 
cies of G and C respectively. In prokaryotic genomes the 
GC skew tends to have a positive value on the leading 
strand of DNA synthesis and a negative value on the 
lagging strand, resulting in polarity changes at the origin 
and terminus of replication (Bentley and Parkhill, 2004). 
The putative position of the replicative terminus is there- 
fore operationally defined as the peak of the cumulative 
GC skew and it typically resides opposite to the origin of 
replication in bacterial genomes (Fig. 2) (Bentley and 
Parkhill, 2004). For T1E, the peak GC skew was indeed 
mapped opposite the replication origin (at 49.2%) of the 
genome. We also found that the oriC site locates between 
the rpmH and dnaA genes and contains two identical 
boxes (5'-TTATCCACA-3') with the first T of the first box 
corresponding to position 991824 while the last A of 
second box is 991893. Gene names were taken from Best 
Blast Hit when available, and gene products were classi- 
fied into COG category (Tatusov etal., 2003), Pfam, Prk 
and Smart families with RPSBIast. Putative ribosomal 
binding sites and tRNA genes were identified with Rfam 
(Griffiths-Jones etal., 2003) and tRNAscan-SE (Lowe 
and Eddy, 1997). Manual validation and visualization of 
the entire metabolic potential of DOT-T1 E was performed 
using the Pathway Tools Program v. 16.0 (http:// 
bioinformatics.ai.sri.com/ptools/) (Letunic etal., 2008; 
Karp etal., 2010), which allows graphic visualization of 
the P. putida annotations. Analyses were performed using 
an Intel(R) Core (TM)i 7-2600 CPU 3.40 GHz with 8 Gb of 
RAM memory running a linux Ubuntu 11.04 operating 
system. Gene products were analysed, compared and 
assigned to metabolic pathways according to the 
MetaCyc scheme (Caspi etal., 2008), and published 
research articles. The cut-off criteria for identifying 
orthologous proteins were compiled by protein-protein 
pairwise analysis and reciprocal tBLASTN analysis to 
identify mutual best hits as potential orthologues. The 
functional annotations of DOT-T1 E genes were corrected 
for consistency with their counterparts in P. putida KT2440 
and P. putida F1 . The coordinates of numerous genes 
were adjusted according to the criteria of full-length align- 
ment, plausible ribosome binding sites, and minimal 
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Fig. 1. Circular genome of Pseudomonas putida DOT-T1E. G+C content and the three tetranucleotide parameters are plotted on the inner- 
most four rings. Distance (second innermost circle) is the distance between global and local sliding window tetranucleotide patterns. Pattern 
skew (third inner most circle) is the distance between tetranucleotide rankings on direct and reverse strands. Oligonucleotide variance (fourth 
inner most circle) is the numerical variance of oligomers, where a lower value indicates tetramer usage and is more highly restricted (for 
example in repeat regions) (Klockgether era/., 2011). The third and second outermost circles show the frequency of distribution of overrepre- 
sented (% 2 > 3000) and highly overrepresented (x 2 > 7000) 8-1 4 mers in the genome of P. putida DOT-T1 E. The outermost ring visualizes dif- 
ferences between tetranucleotide usage and the frequency of the overrepresented longer oligomers. Figures were created with JcircleGraph 
(Davenport era/., 2009). 



overlap between genes on opposite DNA strands. 
Figure 1 shows the Genome Atlas of P. putida (Ussery 
etal., 2009). 

We analysed the genome to identify potential genomic 
islands using three different algorithms based on: (i) lack 
of continuity in the genome, (ii) alignment to other 
P. putida strains and (iii) G+C content and codon usage. 
This yielded four island regions, 1 504 914-1 553 486; 

3 046 659-3 066 609; 4 526 081-4 539 056 and 

4 945 609-4 985 959. Most ORFs in these four islands 
encode hypothetical proteins of unknown function. ORFs 
in islands 1 and 4 exhibit no homology with any other 
known sequence, although significant homology was 
found with transposases. ORFs in island 3 and 2 are 
conserved in P. putida ND6, a strain that degrades naph- 
thalene (Li era/., 2012). 



Metabolic potential 

As indicated above analysis of the entire metabolic poten- 
tial of DOT-T1 E was performed using the Pathway Tools 
Program v.16.0 (http://bioinformatics.ai.sri.com/ptools/) 
(Karp etal., 2002; Letunic era/., 2008). In the genome of 
P. putida DOT-T1E we identified up to 1751 enzymatic 
reactions performed by approximately 1686 enzymes with 
1268 unique potential substrates. A numerical classifica- 
tion for the enzymes based on the chemical reactions they 
carried out according to the Enzyme Commission number 
(EC number) was elaborated in order to understand the 
metabolic potential of this strain. According to EC nomen- 
clature (Bairoch, 2000), oxidoreductases (EC 1) were 
the most abundant enzymes, representing 41% of the 
total (Fig. 3A). Enzymes belonging to EC classes 2 
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Fig. 2. Pseudomonas putida DOT-T1E chromosome GC Skew analysis. Gen Skew is defined as the normalized excess of G over C in a 
given sequence. It is given by (G-C)/(G+C), and it is calculated with a sliding window of 1000 nucleotides along the genome. It is represented 
in blue. The cumulative GC-skew is the sum of the values of neighbouring sliding windows from an arbitrary start to a given point in the 
sequence and it is represented in red. GC-skew is positive in the leading strand and negative in the lagging strand. 



(transferases), EC classes 3 (hydrolases) and 4 (lyases) 
represented 21%, 17% and 10% of all enzymes respec- 
tively, while isomerases (EC 5) and ligases (EC 6) were 
the least abundant, with 5% and 6% of total enzymes 
respectively. This is consistent with the scenario of a 
high metabolic versatility described for Pseudomonads 
(Daniels etal., 2010; Palleroni, 2010). 

The second level of EC nomenclature (EC X.X) 
includes a total of 65 subclasses, of which 51 are present 
in P. putida DOT-T1E (Fig. 3B). As expected, from the 
high number of oxidoreductases, two subclasses of this 
group were among the most abundant with enzymes that 
use the CH-OH group as donor (EC 1.1) and those using 
aldehyde as donors (EC 1.2) representing nearly 12% of 
the total for each group. A striking observation was the 
presence of certain abundant enzyme classes, such as 
for example phosphotransferases (EC 2.7, 7% of total); 
and a series of hydrolases acting on carbon-nitrogen 
bonds (EC 3.5, 5% of total), or acting on ester bonds 
and anhydrides (EC 3.1; about 5% of total). Figure 3B 
presents the enzymes of DOT-T1E grouped based on 
their subclasses. We further classified the enzymes iden- 
tified in functional subclasses according to the EC X.X.X 
nomenclature to focus on the potential donors and accep- 
tors in the case of oxidoreductase enzymes or potential 



groups of substrates in other enzymes (Fig. 3C). Among a 
total number of 269 subclasses in the third level of EC 
nomenclature (EC X.X.X), 150 were present in P. putida 
DOT-T1E. Oxidoreductases using aldehydes as donor 
groups with NAD + or NADP + as acceptor (EC 1 .2.1) were 
the most abundant (11% of the total), also numerically 
important were the carbon-oxygen lyases (EC 4.2.1, 4% 
of total), nucleotidyl phosphotransferases (EC 2.7.7, 3% 
of total) and acyltransferases (EC 2.3.1, 3% of total). 

The enzyme data sets were additionally used to 
analyse potential substrates and to generate a complete 
list of enzyme distribution per functional category EC 
X.X. X.X, the data for which is shown in Table S1. Using 
the Pathway Tool platform, the set of phenomics assays 
previously described by our group (Daniels etal., 2010), 
and the EC X.X.X classification allowed us to explain the 
pattern of growth of strain DOT-T1E with 65 different 
carbon sources, 60 nitrogen sources, and 15 sulfur 
sources used as nutrients (Table S2). In total 425 path- 
ways for metabolism of different compounds were delin- 
eated. This analysis confirms the limited ability of P. putida 
to use sugars as a C source, which is restricted to 
glucose, gluconate and fructose. DOT-T1E has a com- 
plete Entner-Doudoroff route for utilization of glucose and 
other hexoses, but lacks the 6-phosphofructokinase of the 
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Fig. 3. Distribution of enzyme activities of P. putida DOT-T1E classified according to the EC nomenclature. (A) EC X; (B) EC XX; and (C) EC 
XXX. Colour code for classes and subclasses by numbers are indicated. For full details of the EC classification the reader is referred to http:// 
www.chem.qmul.ac.uk/iubmb/enzyme/. 



glycolytic pathway, in agreement with the genome analy- 
sis of others Pseudomonads (del Castillo etal., 2007). A 
large number of sugars were found to not be metabolized 
by T1E including xylulose, xylose, ribulose, lyxose, 
mannose, sorbose, D-mannose, alginate, rhamnose, 
rhamnofuranose, galactose, lactose, epimelibiose, raffi- 
nose, sucrose, stachyose, manninotriose, melibiose, 
tagatose, starch and cello-oligosaccharides, to cite some, 
in agreement with the lack of genes for the metabolism of 
these chemicals after the genome analysis of this strain. 
The results also confirmed the ability of P. putida to use as 
a C source organic acids (such as acetic, citric, glutaric, 
quinic, lactic and succinic among others), certain /.-amino 
acids (Ala, Arg, Asn, Glu, His, He, Lys, Pro, Tyr and Val), 



and various amino organic compounds. (See Figs S1-S4 
for examples of catabolic pathways for sugars, 
amino acids, organic acids and aromatic compounds 
catabolism.) 

Strain T1E harbours genes for a limited number of 
central pathways for metabolism of aromatic compounds 
and numerous peripheral pathways for funnelling of 
aromatic compounds to these central pathways. As in 
other Pseudomonads one of the strategies exploited by 
this microbe for the degradation of different aromatic com- 
pounds is to modify their diverse structures to common 
dihydroxylated intermediates (Dagley, 1971); another 
strategy is to generate acyl-CoA derivatives such as 
phenylacetyl-CoA (Fernandez era/., 2006). Regarding 
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peripheral pathways the P. putida DOT-T1 E genome 
analysis has revealed determinants for putative enzymes 
able to transform a variety of aromatic compounds. The 
DOT-T1E strain is able to use aromatic hydrocarbons 
such as toluene, ethylbenzene, benzene and propylben- 
zene to cite some (Mosqueda era/., 1999). The strain 
also uses aromatic alcohols such as conyferyl- and 
coumaryl-alcohols and their aldehydes; a range of aro- 
matic acids such as ferulate, vanillate, p-coumarate, 
p-hydroxybenzoate, p-hydroxyphenylpyruvate, pheny- 
Ipyruvate, salicylate, gallate and benzoate (see Fig. S4). 
These chemicals are channelled to central catabolic path- 
ways. Upon oxidation of these chemicals they are 
metabolized through one of the three central pathways for 
dihydroxylated aromatic compounds present in this strain. 
The p-ketoadipate pathway is a convergent pathway for 
aromatic compound degradation widely distributed in soil 
bacteria. This pathway consists of a catechol branch {cat) 
and protocatechuate branch (pea). The pea genes in 
P. putida DOT-T1 E are arranged in three operons [pca- 
RKFTBDC (T1E_0230 through T1E_0238), pcaGH 
(T1E_0829 and T1E_830), pcaJI (T1E_2058 and 
T1E_2059)], as is also the case in other P. putida and 
P. syringae strains (Fig. S5). 

The car genes encode the proteins responsible for cat- 
echol degradation and are organized in two clusters 
[catRBCA (T1E_5502 through T1E_5505) and catBCA 
(T1E_1744 through T1E_ 1746)] (Fig. S6), maintaining 
the gene order found in others P. putida strains and also in 
P. aeruginosa. The identity of the catBC and A genes in 
both clusters is in the range of 79-82%. In addition, we 
should mention that two other catA genes were found, 
one of them with a high degree of similarity to the KT2440 
catA2 gene, which corresponded to ORF T1 E_1 057, that 
is adjacent to the benRABCDK genes (T1E_1055 to 
T1 E_1 064) for benzoate degradation; while the other catA 
allele corresponded to ORF T1 E_551 1 . It should be noted 
that this allele is within a cluster of genes that are tran- 
scribed in the same direction and which encode genes for 
salycilate metabolism (T1E_5510 through T1E_5513). 

The genes involved in phenylacetate degradation were 
also identified in P. putida DOT-T1E. There are 16 genes 
encoding for phenylacetate degradation organized in a 
cluster (ORFs T1E_5587 to T1E_5603) and within the 
cluster a series of potential operons were identified, i.e. 
the paaGHIJK genes (T1 E_5590 through T1 E_5594) that 
encode the ring-hydroxylating oxygenase enzyme, the 
paaABCDE genes that encode the (3-oxidation enzymes, 
a potential phenylacetate transport system (paaLM) and 
the regulatory system made of paaXY, that correspond to 
T1 E_5587 and T1 E_5588 respectively. 

Homologous genes for degradation of homogentisate 
are also present in strain DOT-T1E. Homogentisate is 
catabolized by a central catabolic pathway that involves 



t_i M A llu 6.3.4.6 (_, M A..,J^ ftu 3.5.1.54 m 

H2N NH2 H2N N OH »- h H 

.,„_ TIE 3118 , , . TIE 3119 

urea tie~3809 urea-l-carboxylate T1EI3808 ammonia 



tie>3°5 ammonia 

T1E_4306 

Fig. 4. Pathway for utilization of urea as an N source by P. putida. 
The genes that encoded the enzymes of these two pathways were 
identified based on blast analysis and comparison to proteins that 
carry out the indicated reactions. 

three enzymes, homogentisate dioxygenase (T1E_1557), 
a newly identified putative maleylacetoacetate isomer- 
ase (T1E_1555) and fumarylacetoacetate hydrolase 
(T1E_1558). In this pathway homogentisate is funnelled 
to yield fumarate and acetoacetate. A search for hpa and 
gtd genes that encode genes belonging to the homopro- 
tocatechuate and gentisate pathways yielded no results 
from the DOT-T1 E genome, which suggests the absence 
of a meta ring-cleavage pathway for the degradation of 
homoprotocatechuate and gentisate. 

Pseudomonads strains are able to use a range of inor- 
ganic nitrogen sources. In this regard three predicted 
transporters involved in the uptake of ammonium were 
identified. T1E incorporates ammonium into C skeletons 
using mainly the ATP-dependent activity of glutamine syn- 
thetase (GS) followed by the action of glutamate synthase 
(GOGAT). The genome of T1E encodes four GS 
(T1E_0118, 1260, 2050 and 4444) and four GOGAT 
enzymes (T1E_1644, 2053, 2506 and 3293). Strain T1E 
can use nitrate as an N source, which is reduced to 
ammonium using an assimilatory nitrate reductase (EC: 
1.7.99.4) encoded by the T1E_4793 gene, that is in a 
cluster with nirB and nirD which encode an assimilatory 
nitrite reductase (EC1.7.1.4). (The ORFs encoding these 
proteins correspond to T1E_4793 through T1E_4795.) 
The strain also has the complement of genes for utiliza- 
tion of urea either through direct conversion to ammonia 
(T1E_4304 through T1E_4306, ureABC) or via conver- 
sion first to urea-1-carboxylate (T1E_3118 through and 
3809) and then conversion to ammonia (T1E_3119 and 
T1E_3808) (Fig. 4). 

Details for the utilization of D- and /.-amino acids as N 
sources were published by Daniels and colleagues 
(2010). It was found that the wild-type DOT-T1E strain 
was able to use a number of either D- or /.-amino acids 
(i.e. D-ornithine, D-alanine, D-arginine, D-asparagine, 
D-lysine and D-valine), dipeptides, ethanolamine, and 
adenine as an N source (Daniels era/., 2010). It is of 
interest to highlight that this strain can use several 
D-amino acids for which racemases are needed. We have 
found that the genome of DOT-T1E encodes at least 
five broad-substrate racemases (T1E_2780, TIE_3429, 
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TIE_1731, TIE_0166, TIE_4880) that can convert 
D-amino acids into /.-amino acids which upon transami- 
nation allow the catabolism of these compounds to 
provide nitrogen for growth (Daniels era/., 2010). Eight 
aminopeptidases (TIE_3567, TIE_2564, TIE_4792, 
TIE_1957, TIE_2243, TIE_3241, TIE_3898, TIE_0833) 
also allows this bacterium to utilize a number of dipeptides 
and tripeptides as C- and N- sources, in agreement with 
the saprophytic character of strains of this species 
(Daniels etal., 2010). 

Strain T1E has a number of genes that may encode 
enzymes/transporters needed for the acquisition of inor- 
ganic phosphate, namely: (i) two low-affinity Pit type 
transporters (T1 E_0227 and T1 E_0045), (ii) two putative 
ABC-type inorganic phosphate high-affinity transporter 
(T1E_2661 through 2663 and T1E_3987 through 3989) 
and (iii) a PstS type (T1 E_2660) high-affinity transporter 
system regulated by the phoBR (T1E_3994 and 3993) 
response regulator system. This strain uses organic phos- 
phate ester compounds under phosphorous-limiting 
conditions (Daniels etal., 2010). T1E also use organic 
phosphonates that are transported by a high-affinity ABC 
transport system consisting of the phnD, phnE and phnC 
gene products (T1E_4609 through 4612). 

Members of the pseudomonadaceae have been 
reported to play a key role in mineralization of carbon 
bound sulfur in rhizosphere soils. Organic sulfur in soils is 
comprised mostly of sulfonates and sulfate esters; hence, 
many soil bacteria carry genes that encode enzymes for 
utilization of alkanesulfonates. Metabolism of these com- 
pounds is achieved through the action of the Ssu enzymes, 
which are encoded by a set of genes that form an operon, 
namely, ssuA through F(T1E_2976 through 2982). This 
organization is similar to that in other Pseudomonas 
(Kahnert and Kertesz, 2000). The strain DOT-T1E is also 
endowed with at least one putative arylsulfatase 
(T1 E_5507) which may explain the ability of the strain to 
use aromatic sulfate esters (Daniels etal., 2010). The 
DOT-T1 E strain is also endowed with four genes that may 
encode the enzymes required to make sulfur available 
from methionine (T1 E_0568, T1 E_2981 , T1 E_4829 and 
T1 E_4830), which is released as sulfite (Fig. S7). The set 
of reactions is initiated by MdeAas in other pseudomonads 
and the pathway is depicted in Fig. S7. 

A relevant characteristic of DOT-T1 E is its capability to 
grow on minimal medium without the need of vitamins or 
other cofactors. We found 165 genes encoding enzymatic 
reactions mediating the biosynthesis of a number of 
cofactors, i.e. nicotinate, nicotinamide, vitamin B6, ribofla- 
vin, ubiquinone, porphyrin, biotin, thiamine, folate, pan- 
tothenate and CoA which amounts for 74 distinct 
biosynthetic pathways. This is consistent with a metabo- 
lism in which different enzymes have been described to 
use these molecules as cofactors. 



Based on phenotypic analysis using the BIOSCREEN 
growth test system described by Daniels and colleagues 
(2010), it was shown that P. putida T1 E tolerated various 
heavy metals. Based on the strain's genome sequence, 
64 genes were identified that encode proteins putatively 
involved in heavy metal resistance and homeostasis 
(Table 1). The majority of the P. putida T1E heavy metal 
resistance genes are found spread throughout the 
genome, and they are conserved among all sequenced 
P. putida strains. 

Up to three different systems potentially involved in 
simultaneous cobalt, zinc and cadmium resistance were 
found. One of the cation efflux systems is the CzcD 
(T1 E_2808) immersed in a cluster with the corresponding 
response regulator CzcR (T1E_2811) and the sensor his- 
tidine kinase encoded by the czcS gene (T1E_2812). 
Another family of transporters that may mediate the extru- 
sion of these three heavy metal ions are the one encoded 
by the cadA 1 (T1 E_2820) and cadA2 (T1 E_4489) genes; 
as well as by the resistance-nodulation-cell division 
(RND) pump CzcABC (T1E_5270, T1E_5271, T1E_ 
5272). The CusABC efflux system (T1E_4694, T1E_ 
4695, T1E_4696) is involved resistance to silver and 
copper ions. Seven genes involved in resistance to 
arsenite-arsenate-antimonite efflux were annotated. 
Four of them arsHCBR made an operon (T1E_2719- 
2722), and the three other genes related to arsenite 
resistance (T1E_4939, T1E_4996 and T1E_1144) are 
scattered throughout the genome. Finally one chromate 
resistance protein ChrA (T1E_3354) was found in the 
genome of T1 E suggesting it is the responsible for chro- 
mate efflux in this strain. 

Biotransformation potential 

As mentioned above DOT-T1 E has the ability to thrive in 
the presence of toxic organic solvents that normally form 
a biphasic system with water. This property can be 
exploited to develop double-phase biotransformation 
systems (organic solvent and water) in which water 
insoluble chemicals, toxic substrates or chemical prod- 
ucts are kept in the organic phase. The main advantages 
of these systems are that the product(s) is(are) continu- 
ously removed by a solvent phase, their toxic effects are 
decreased and the lifespan of the biocatalytic system is 
longer. In addition, if the concentration of the product 
increases in the organic phase, product recovery is easier 
and less costly (Bruce and Daugulis, 1991; Leon etal., 
1998). Rojas and colleagues (2004) demonstrated that 
P. putida DOT-T1 E was tolerant to different aliphatic alco- 
hols such as decanol, nonanol and octanol. These 
aliphatic alcohols are useful in double-phase biotransfor- 
mation systems to deliver hydrophobic or toxic com- 
pounds or to recover added value products that partition 
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Table 1. Proteins found in P. putida DOT-T1E that are associated with metal resistance and homeostasis. 



Gene Protein 



location 


name 


Metal 


Family/domain 


Predicted role 


Definition 


E-value 


T1E 


_0296 


CzcS2 


Me 2 


TC reg 


Sensor protein 


Hypothetical protein 


0.0 


T1E_ 


_0297 


CzcR2 


Me 2 


TC reg 


Response regulator 


DNA-binding response regulator CzrR 


e-126 


T1E. 


_0503 


TPMT 


Te, Se 


TPMT 


Te and Se Se/Te 
detoxification 


Thiopurine S-methyltransferase 


e-122 


T1E. 


_0621 


ZnuA2 


Zn/Mn(?) 


PBD 


Zn/Mn(?) uptake 


Periplasmic solute-binding protein 


e-167 


T1E_ 


_0622 


ZnuC2 


Zn/Mn(?) 


ATP-binding protein 


Zn/Mn(?) uptake 


Cation ABC transporter, AP-binding protein 


e-122 


T1E_ 


_0658 


Fur 


Fe 


Fur 


Fe regulation 


Ferric uptake regulator, Fur family 


7e-073 


T1E_ 


_0727 


CopA 


Cu 


HMA 


Copper exporting ATPase 


Heavy metal translocating P-type ATPase 


0.0 


T1E 


_1144 


ArsR3 


As, Sb 


ArsR 


Transcriptional regulator 


ArsR family transcriptional regulator 


5e-152 


T1E_ 


J 232 


CopS2 


Cu 


TC reg 


Response regulator 


Heavy metal sensor signal transduction histidine kinase 


0.0 


T1E, 


_1 233 


CopR2 


Cu 


TC reg 


Sensor protein 


Two-component heavy metal response transcriptional 
regulator 


e-122 


T1E. 


J 234 


T1E_1234 


Cu 


Cupredoxine 


Copper homeostasis 


Plastocyanin/azurin family copper-binding protein 


4e-072 


T1E. 


_1474 


ModC 


Mo 


ATP-binding protein 


Mo uptake 


Molybdate ABC transporter ATPase 


0.0 


T1E_ 


J 475 


ModB 


Mo 


I M pore 


Mo uptake 


Molybdate ABC transporter inner membrane protein 


e-128 


T1E, 


J 476 


ModA 


Mo 


PBD 


Mo uptake 


Molybdenum ABC transporter periplasmic 
molybdate-binding protein 


e-136 


T1E 


J 824 


NikE 


Ni 


ATP-binding protein 


Ni uptake 


Nickel transporter ATP-binding protein NikE 


3e-041 


T1E_ 


J 827 


NikB 


Ni 


I M pore 


Ni uptake 


Nickel transporter permease NikB 


7e-050 


T1E_ 


_2011 


CumA 


Cu 


Cu oxidase 


Copper homeostasis 


Multicopper oxidase 


0.0 


T1E, 


_2070 


NikC 


Ni 


I M pore 


Ni uptake 


Nickel transporter permease NikC 


4e-059 


T1E_ 


_2193 


ModR 


Mo 


modE 


Mo uptake regulation 


Mode family transcriptional regulator 


e-127 


T1E_ 


_2274 


T1 E_2274 


Cu 


(MFS) transporters 


Copper homeostasis 


Bcr/CfIA family multidrug resistance transporter 


0.0 


T1E, 


_2279 


T1 E 2279 


Co/Zn/Cd 


OEP 


Cobalt-zinc-cadmium 
resistance 


Heavy metal RND efflux outer membrane protein, CzcC 
family 


0.0 


T1E 


_2577 


ZnuB2 


Zn/Mn(?) 


I M pore 


Zn/Mn uptake 


hypothetical protein 


e-160 


T1E_ 


_2719 


ArsH2 


As, Sb 


ArsH 


Arsenical resistant 


ArsH protein 


e-132 


T1E 


_2720 


ArsC2 


As, Sb 


ArsC 


As(V) reduction 


Arsenate reductase 


8e-086 


T1E_ 


_2721 


ArsB2 


As, Sb 


ArsB 


As(lll), Sb(lll) efflux 


Arsenite efflux transporter 


0.0 


T1E_ 


_2722 


ArsR2 


As, Sb 


ArsR 


Transcriptional repressor 


Arsenic resistance transcriptional regulator 


8e-063 


T1E_ 


_2794 


NikA 


Ni 


PBR 


Ni uptake 


Nickel ABC transporter, periplasmic nickel-binding protein 


1 e-039 


T1E. 


_2808 


CzcD 


Me 2 


CDF 


Transport and regulation 


CDF family cobalt/cadmium/zinc transporter 


e-166 


T1E_ 


_2811 


CzcR1 


Me 2 


TC reg 


Response regulator 


DNA-binding response regulator CzrR 


e-124 


T1E_ 


_2812 


CzcS1 


Me 2 


TC reg 


Sensor kinase 


Sensor histidine kinase 


0.0 


T1E_ 


_2820 


CadA1 


Zn/Cd 


P-type ATPase 


Me2+ efflux 


Heavy metal translocating P-type ATPase 


0.0 


T1E 


_2933 


TetR 


Drug(?) 


TetR 


Transcriptional regulator 


TetR family transcriptional regulator 


e-117 


T1E_ 


J3354 


ChrA 


Cr 


ChrA 


Chromate efflux 


Chromate transporter 


0.0 


T1E_ 


_3756 


PacR(CueR) 


Cu/Ag 


MerR 


Transcriptional regulaor 


MerR family transcriptional regulaor 


1 e-075 


T1E 


_3757 


PacS 


Cu 


P-type 


Cooper uptake 


Heavy metal translocating P-type ATPase 


0.0 


T1E. 


_3759 


PacZ(CopZ) 


Cu 


HMA 


Activator 


heavy metal transport/detoxification protein 


7e-031 


T1E. 


_3760 




Cu 




Cooper homeostasis 


Multidrug resistance transporter, Bcr/CfIA family 


1e-163 


T1E_ 


_4452 


Mfpll 


Me 2 /drug 


RND MFP/HlyD 


Me2+/drug efflux 


efflux transporter, RND family, MFP subunit 


0.0 


T1E_ 


_4453 


Mfpl 


Me 2 /drug 


RND MFP/HlyD 


Me2/drug efflux 


RND efflux transporter 


0.0 


T1E 


_4454 


CzcA4 


Me 2 /drug 


RND 


Me2/drug efflux 


Acriflavin resistance protein 


0.0 


T1E_ 


_4488 


CadR 


Zn/Cd 


MerR 


Cd, Zn efflux 


MerR family transcriptional regulator 


8e-082 


T1E_ 


_4489 


CadA2 


Cd/Zn 


P-type ATPase 


Cd, Zn efflux 


Heavy metal translocating P-type ATPase 


0.0 


T1E_ 


_4513 


CopA1 


Cu 


MultiCU oxidases 


Cu chelation 


Copper resistance protein A 


0.0 


T1E, 


_4672 


ZnuC1 


Zn 


ATP-binding protein 


Zn uptake 


Zinc ABC transporter ATP-binding protein 


e-147 


T1E, 


_4694 


CusA 


Me 2 


RND 


Me2 efflux 


CzcA family cobalt/zinc/cadmium efflux transporter 
permease 


0.0 


T1E. 


_4695 


CusB 


Me 2 


RND MFP/HlyD 


Me2+ efflux 


CzcB family cobalt/zinc/cadmium efflux transporter 

momhrano fiiQir^n nrAtoin 

IIIUIIILHdllt? 1 UolUI 1 LflULcN 1 


0.0 


T1E. 


_4696 


CusC 


Me 2 


OEP 


Me2+ efflux 


CzcC family cobalt/zinc/cadmium efflux transporter outer 
membrane protein 


0.0 


T1E 


_4697 


PorD 




Porin 


Channel basic amino acids 


Porin, putative 


0.0 


T1E_ 


_4698 


CzcR3 


Me 2 


TC reg 


Response regulator 


DNA-binding heavy metal response regulator putative 


e-127 


T1E_ 


_4760 


ZnuA1 


Zn 


PBD 


Zn uptake 


Periplasmic solute-binding protein 


e-171 


T1E_ 


_4761 


Zur 


Zn 


Fur 


Regulator 


FUR family transcriptional regulator 


2e-074 


T1E. 


_4763 


ZnuB1 


Zn 


I M pore 


Zn uptake 


Hypothetical protein 


e-138 


T1E_ 


_4936 




Cu 


CBS 


Copper homeostasis 


CBS domain containing protein 


1e-130 


T1E_ 


_4939 


ArsR1 


As, Sb 


ArsR 


Transcriptional regulator 


ArsR family transcriptional regulator 


5e-048 


T1E_ 


_4996 


ArsC3 


As, Sb 




Arsenate reductase 


Arsenate reductase 


1e-060 


T1E 


_5088 


CzcA5 


Me 2 /drug 
Me 2 /drug 


RND 


Me2+/drug efflux 


Acriflavin resistance protein 


0.0 


T1E_ 


_5089 


Mtrc2 


RND MFP/HlyD 


Me2+/drug efflux 


Efflux transporter, RD family, MFP subunit 


0.0 


T1E_ 


_5270 


CzcA2 


Me 2 


RND 


Cation efflux 


Cobalt-zinc-cadmium resistance protein CzcA 


0.0 


T1E 


_5271 


CzcB2 


Me 2 


RND MFP/HlyD 


Cation efflux 


Cobalt-zinc-cadmium resistance protein CzcB, putative 


0.0 


T1E. 


_5272 


CzcC2 


Me 2 


OEP 


Cation efflux 


Cobalt-zinc-cadmium resistance protein CzcC, putative 


0.0 


T1E. 


_5277 


ompR 




Response regulator 


Copper homeostasis, 
Cobalt-zinc-cadmium 
resistance 


DNA-binding heavy metal response regulator 


6e-124 


T1E. 


_5753 


CopB1 


Cu 


OM protein 


Cu chelation(?) 


Copper resistance B precursor 


e-159 
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preferentially in the organic phase. This concept was 
exploited by Neumann and colleagues (2005) who 
showed that DOT-T1E in the presence of 1-decanol tol- 
erated up to 200-fold higher concentrations of the model 
substrate 3-nitrotoluene than in aqueous medium. In the 
same line Wierckx and colleagues (2008) showed that 
phenol production from glucose by P. putida S12, another 
solvent tolerant strain, increased up to 10-fold using a 
biphasic system. This set of results is the basis that 
support the potential of DOT-T1 E as a useful biocatalyst 
for biphasic systems. 

As described above, a wide range of oxido-reductase 
enzymes are encoded in the genome of DOT-T1E, a 
number of which are of commercial interest. Among these 
are a numerous dioxygenases that might selectively 
hydroxylate the aromatic rings at positions 1 and 2; 2 and 
3; 2 and 5; 3 and 4; and 4 and 5 (see Table S3). These 
dioxygenases may catalyse the stereo-specific dioxy- 
genation of hydrocarbons and could yield secondary com- 
modity chemicals such as adipic acid and y-caprolactam. 

At least 16 monooxygenases that may act on diverse 
chemicals have also been annotated (Table S4). Some of 
these enzymes have the potential to oxidize alkanes to 
their pertinent alcohols, and are of interest to generate 
added-value products such as linear branched alcohols, 
aromatic alcohols, diols, hydroxypropionic acid and 
others. Since the produced chemicals are not metabo- 
lized by P. putida DOT-T1E, they accumulate in culture 
supernatants and high yields can be achieved via extrac- 
tion of the second phase (Fig. S8). 

Biotransformations based on genetically engineering 
production strains for certain compounds may require 
either blocking an existing pathway, recruiting new 
enzymes, or a combination of both approaches. Some of 
these approaches have been used before with DOT-T1 E 
or related solvent-tolerant strains. An example of biotrans- 
formation through inhibition of a single gene is the pro- 
duction of 3-methylcatechol from toluene. Thus, Rojas 
and colleagues (2004) showed that a catechol 2,3- 
dioxygenase knockout mutant of T1 E in a 1-octanol/water 
bioreactor produced 20-fold higher amounts of the com- 
pound than in the aqueous medium; these results dem- 
onstrate the usefulness of double-phase systems. 

Ramos-Gonzalez and colleagues (2001) developed 
a system for transformation of toluene into 4- 
hydroxybenzoate which involves the use of a double 
mutant of T1E in which the toluene dioxygenase and 
p-hydroxybenzoate hydroxylase genes were first inacti- 
vated. Then, a set of genes for sequential oxidation of 
toluene via toluene 4-monooxygenase were incorporated 
and the recombinant strain system produced up to 35 mM 
of the product. 

Efficient bioconversion of glucose to phenol or cinnamic 
acid was achieved by Wierckx and colleagues (2005) and 



Nijkamp and colleagues (2005), respectively, with 
P. putida S12 to this end a tyrosine phenol lyase (EC 
4.1 .9.22) from Pantoea agglomerans was recruited and it 
enabled the S12 strain to produce phenol (Fig. 5). In 
fed-batch assays in water, the productivity was limited by 
accumulation of 5 mM phenol in the medium; above this 
concentration phenol was toxic. However, this toxicity was 
overcome by use of 1-octanol as a second phase and as 
an extractant for phenol in a biphasic system. This 
approach resulted in accumulation of nearly 50 mM 
phenol in the octanol phase (Wierckx era/., 2008). Other 
possibilities for the production of added-value molecules 
with DOT-T1E are their synthesis from tyrosine; for 
example, DOT-T1 E can produce L-DOPA from tyrosine 
(Fig. 5). This can be achieved by recruiting one of the 
following activities: a polyphenol oxidase (EC 1.10.3.1), a 
tyrosinase (1 .14.18.1) or a tyrosine 3-monooxygenase (E 
1 .1 4.1 6.2) (Krishnaveni et al., 2009; Surwase and Jadhav, 
2011). It should also be noted that with tyrosine as a 
substrate DOT-T1 E can produce tyramine (via an internal 
aromatic amino acid decarboxylase, EC 4.1.1.28) and 
4-hydroxyphenylpyruvate using a tyrosine amino trans- 
ferase (EC 2.6.1.5). However, accumulation of the prod- 
ucts of these biotransformations requires the inhibition of 
further catabolism of the products because they can be 
used as a C source by DOT-T1E (Daniels era/., 2010). 

One of our aims is to customize strains for the produc- 
tion of aromatic alcohols for biofuel production. In this 
regard DOT-T1E can be used to produce alkyl and aro- 
matic alcohols (Fig. 6) through blocking the catabolic 
pathways for amino acid degradation, in which keto acid 
intermediates are converted into their corresponding alco- 
hols, as reported for Escherichia coli (Atsumi era/., 2008) 
- a process that requires the recruitment of a keto acid 
decarboxylase to produce an intermediate aldehyde 
that is subsequently transformed into its corresponding 
alcohol. 

D-xylose is the second most abundant sugar in ligno- 
cellulosic materials and its utilization by industrial organ- 
isms to produce biofuels and added-value aromatic 
compounds is of interest (Octave and Thomas, 2009). As 
described above, strains of the P. putida species cannot 
use pentose sugars, but this was overcome via the engi- 
neered addition of xylAB genes, which allow the conver- 
sion of D-xylose in D-xylulose and xylulose-5-P, to allow 
metabolism of D-xylose via the pentose phosphate 
pathway (Meijnen et al., 2008; 2009). However, growth 
was rather slow and fast growers were isolated after 
enrichment in fermentors. In a recent omics-based study, 
the authors have shown that high yield growth involved 
inactivation of glucose dehydrogenase and rearrange- 
ment of central carbon catabolism to allow for more effi- 
cient decarboxylation of 6-phosphogluconate for the 
catabolism of the sugar via the pentose phosphate 
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Homogentisate 
k 




Fig. 5. Biotransformation of tyrosine by 
P. putida through metabolic blockage or gene 
recruitment. The EC XXXX of the enzymes 
needed for the listed biotransformation are 
indicated. The text describes the approaches 
used by different research groups to achieve 
the indicated products. 
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OH 
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recruiting 
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TIE 3359 
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HO 



NH 2 



Tyramine 



T1EJJ35* =CaSa__ 



4-Hydroxyphenyl 
acetaldehyde 



pathway (Meijnen etal., 2012). Since P. putida S12 is as 
tolerant to solvents as DOT-T1 E (Segura et al., 2003), we 
hypothesize that a similar strain of DOT-T1E could be 
engineered. 



In summary, the analysis of the genome of the solvent- 
tolerant DOT-T1 E strain explains the catabolic potential 
of this microorganism in accordance with previously pub- 
lished physiological studies. The use of the Pathway Tool 
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Intermediary 
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2-Ketobutyrate 
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OH 
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Fig. 6. Potential synthesis of different 
alcohols from keto acids by D0T-T1 E. 
2-Ketoisovalerate, 2-ketobutyrate and phe- 
nylpyruvate are produced in the catabolism of 
isoleucine, threonine and tryptophane respec- 
tively. According to Atsumi and colleagues 
(2008) recruitment of a broad substrate range 
keto acid decarboxylase (KDH) yields an 
aldehyde, which along with one of the multiple 
alcohol dehydrogenase enzymes encoded in 
the genome of this strain can lead to the syn- 
thesis of the corresponding alcohol (see 
Table S1). 



o 

Phenylpyruvate Intermediary 2-Phenyl ethanol 
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platform together with the identification of enzymes using 
the international EC codes not only support the metabolic 
reactions, but also provide an opportunity to design 
biotransformation reactions to produce value-added prod- 
ucts in high concentration. Due to the proven ability of 
DOT-T1E to thrive in the presence of a second organic 
phase in biphasic systems, only minimal genetic manipu- 
lation will be required in order to reap substantial reward 
from the genomic analysis reported here. 

Acknowledgements 

We thank Manuel Ferrer for early collaboration in this 
project and M. Mar Fandila for secretarial assistance. 

Conflict of interest 

None declared. 



References 

Arias, S., Olivera, E.R., Arcos, M., Naharro, G., and Luengo, 
J.M. (2008) Genetic analyses and molecular characteriza- 
tion of the pathways involved in the conversion of 
2-phenylethylamine and 2-phenylethanol into phenylacetic 
acid in Pseudomonas putida U. Environ Microbiol 10: 413- 
432. 

Atsumi, S., Cann, A.F., Connor, M.R., Shen, C.R., Smith, 
K.M., Brynildsen, M.P., era/. (2008) Metabolic engineering 
of Escherichia coli for 1-butanol production. Metab Eng 10: 
305-311. 

Bairoch, A. (2000) The ENZYME database. Nucleic Acids 
Res 28: 304-305. 

Bentley, S.D., and Parkhill, J. (2004) Comparative genomic 
structure of prokaryotes. Annu Rev Genet 38: 771-792. 

Bernal, P., Munoz-Rojas, J., Hurtado, A., Ramos, J.L., and 
Segura, A. (2007) A Pseudomonas putida cardiolipin 
synthesis mutant exhibits increased sensitivity to drugs 
related to transport functionality. Environ Microbiol 9: 
1135-1145. 

Bruce, L.J., and Daugulis, A.J. (1991) Solvent selection 
strategies for extractive biocatalysis. Biotechnol Prog 7: 
116-124. 

Caballero, A., Esteve-Nunez, A., Zylstra, G.J., and Ramos, 
J.L. (2005) Assimilation of nitrogen from nitrite and trinitro- 
toluene in Pseudomonas putida JLR11. J Bacteriol 187: 
396-399. 

Caspi, R., Foerster, H., Fulcher, C.A., Kaipa, P., Krummen- 
acker, M., Latendresse, M., era/. (2008) The MetaCyc 
database of metabolic pathways and enzymes and the 
BioCyc collection of pathway/genome databases. Nucleic 
Acids Res 36: D623-D631 . 

Dagley, S.C (1971) Catabolism of aromatic compounds by 
micro-organism. Adv Microb Physiol 6: 1-46. 

Daniels, C, Godoy, P., Duque, E., Molina-Henares, M.A., de 
la Torre, J., del Arco, J.M., era/. (2010) Global regulation of 
food supply by Pseudomonas putida DOT-T1 E. J Bacteriol 
192: 2169-2181. 



Davenport, C.F., Wiehlmann, L, Reva, O.N., and Tiimmler, 
B. (2009) Visualization of Pseudomonas genomic structure 
by abundant 8-14 mer oligonucleotides. Environ Microbiol 
11: 1092-1104. 

Del Castillo, T, Ramos, J.L., Rodriguez-Herva, J. J., Fuhrer, 
T, Sauer, U., and Duque, E. (2007) Convergent peripheral 
pathways catalyze initial glucose catabolism in Pseu- 
domonas putida: genomic and flux analysis. J Bacteriol 
189: 5142-5152. 

Delcher, A.L., Harmon, D., Kasif, S., White, 0., and Salzberg, 
S.L. (1999) Improved microbial gene identification with 
GLIMMER. Nucleic Acids Res 27: 4636-4641. 

Dominguez-Cuevas, P., Gonzalez-Pastor, J.E., Marques, S., 
Ramos, J.L., and de Lorenzo, V. (2006) Transcriptional 
tradeoff between metabolic and stress-response programs 
in Pseudomonas putida KT2440 cells exposed to toluene. 
J Biol Chem 281: 11981-11991. 

Fernandez, C, Ferrandez, A., Mihambres, B., Diaz, E., and 
Garcia, J.L. (2006) Genetic characterization of the 
phenylacetyl-coencime a oxygenase from the aerobic phe- 
nylacetic acid degradation pathway of Escherichia coli. 
Appl Environ Microbiol 72: 7422-7426. 

Garcia, V., Godoy, P., Daniels, C, Hurtado, A., Ramos, J.L., 
and Segura, A. (2010) Functional analysis of new trans- 
porters involved in stress tolerance in Pseudomonas putida 
DOT-T1 E. Environ Microbiol Rep 2: 389-395. 

Godoy, P., Molina-Henares, A.J., Torre, J., Duque, E., and 
Ramos, J.L. (2010) Characterization of the RND family of 
multidrug efflux pumps: in silico to in vivo confirmation of 
four functionally distinct subgroups. Microb Biotechnol 3: 
691-700. 

Griffiths-Jones, S., Bateman, A., Marshall, M., Khanna, A., 
Eddy S.R. (2003) Rfam: an RNA family database. Nucleic 
Acids Res 31: 439-441. 

Heipieper, H.J., de Ward, P., van der Meer, P., Killian, J.A., 
Isken, S., de Bont, J.A.M., et al. (2001 ) Regiospecific effect 
of 1-octanol on cis-trans isomerization of unsaturated fatty 
acids in the solvent-tolerant strain Pseudomonas putida 
S12. Appl Microbiol Biotechnol 57: 541-547. 

Heipieper, H.J., Meinhardt, F, and Segura, A. (2003) The 
cis-trans isomerase of unsaturated fatty acids in Pseu- 
domonas and Vibrio: biochemistry, molecular biology and 
physiological function of a unique stress adaptive mecha- 
nism. FEMS Microbiol Lett 229: 1-7. 

Isken, S., and de Bont, J. A.M. (1996) Active efflux of toluene in 
a solvent-resistant bacterium. J Bacteriol 178: 6056-6058. 

Jimenez, J.L., Mihambres, B., Garcia, J.L., and Di'az, E. 
(2002) Genomic analysis of the aromatic catabolic path- 
ways from Pseudomonas putida KT2440. Environ Micro- 
biol 4: 824-841. 

Junker, F, and Ramos, J.L. (1999) Involvement of the cis/ 
trans isomerase Cti in solvent resistance of Pseudomonas 
putida DOT-T1 E. J Bacteriol 181 : 5693-5700. 

Kahnert, A., and Kertesz, M.A. (2000) Characterization 
of a sulfur-regulated oxygenative alkylsulfatase from 
Pseudomonas putida S-313. J Biol Chem 275: 31661- 
31667. 

Karp, P.D., Paley, S., and Romero, P. (2002) The Pathway 
Tools software. Bioinformatics 18 (Suppl. 1): S225-S232. 

Karp, P.D., Paley, S.M., Krummenacker, M., Latendresse, M., 
Dale, J.M., Lee, T.J., era/. (2010) Pathway Tools version 



© 2013 The Authors. Microbial Biotechnology published by John Wiley & Sons Ltd and Society for Applied Microbiology, Microbial 

Biotechnology 6, 598-611 



Solvent tolerance strategies 609 



13.0: integrated software for pathway/genome informatics 
and systems biology. Brief Bioinform 11: 40-79. 

Keweloh, H., and Heipieper, H.J. (1996) Trans unsaturated 
fatty acids in bacteria. Lipids 31: 129-137. 

Kieboom, J., Dennis, J. J., de Bont, J. A.M., and Zylstra, G.J. 
(1998) Identification and molecular characterization of an 
efflux pump involved in Pseudomonas putida S12 solvent 
tolerance. J Biol Chem 273: 85-91 . 

Kim, K., Lee, S., Lee, K., and Lim, D. (1998) Isolation and 
characterization of toluene-sensitive mutants from the 
toluene-resistant bacterium Pseudomonas putida GM73. 
J Bacteriol 180: 3692-3696. 

Klockgether, J., Cramer, N., Wiehlmann, L., Davenport, C.F., 
and Tummler, B. (2011) Pseudomonas aeruginosa 
genomic structure and diversity. Front Microbiol!: 150. 

Krishnaveni, R., Rathod, V., Thakur, M.S., and Neelgund, Y.F. 
(2009) Transformation of /.-tyrosine to L-DOPA by a novel 
fungus, Acremonium rutilum, under submerged fermenta- 
tion. Curr Microbiol 58: 122-128. 

Leon, R., Fernandes, P., Pinheiro, H.M., and Cabral, J. M.S. 

(1998) Whole-cell biocatalysis in organic media. Enzyme 
Microb Technol 23: 483-500. 

Letunic, I., Yamada, T., Kanehisa, M., and Bork, P. (2008) 
iPath: interactive exploration of biochemical pathways and 
networks. Trends Biochem Sci 33: 101-103. 

Li, S., Zhao, H., Li, Y., Niu, S., and Cai, B. (2012) Complete 
genome sequence of the naphthalene-degrading Pseu- 
domonas putida strain ND6. J Bacteriol 194: 5154-5155. 

Lowe, T.M., and Eddy, S.R. (1997) tRNAscan-SE: a program 
for improved detection of transfer RNA genes in genomic 
sequence. Nucleic Acids Res 25: 955-964. 

Meijnen, J. P., de Winde, J.H., and Ruijssenaars, H.J. (2008) 
Engineering Pseudomonas putida S12 for efficient utiliza- 
tion of D-xylose and L-arabinose. Appl Environ Microbiol 
74: 5031-5037. 

Meijnen, J. P., de Winde, J.H., and Ruijssenaars, H.J. (2009) 
Establishment of oxidative D-xylose metabolism in Pseu- 
domonas putida S12. Appl Environ Microbiol 75: 2784- 
2791. 

Meijnen, J. -P., and de Winde, J.H., and Ruijssenaars, H.J. 
(2012) Metabolic and regulatory rearrangements underly- 
ing efficient D-xylose utilization in engineered Pseu- 
domonas putida S12. J Biol Chem 287: 14608-14614. 

Molina, L., Duque, E., Gomez, M.J., Krell, T., Lacal, J., 
Garcfa-Puente, A., etal. (2011) The pGRT1 plasmid of 
Pseudomonas putida DOT-T1 E encodes functions relevant 
for survival under harsh conditions in the environment. 
Environ Microbiol 13: 2315-2327. 

Mosqueda, G., and Ramos, J.L. (2000) A set of genes encod- 
ing a second toluene efflux system in Pseudomonas putida 
DOT-T1 E is linked to the todgenes for toluene metabolism. 
J Bacteriol 182: 937-943. 

Mosqueda, G., Ramos-Gonzalez, M.I., and Ramos, J.L. 

(1999) Toluene metabolism by the solvent-tolerant Pseu- 
domonas putida DOT-T1 strain, and its role in solvent 
impermeabilization. Gene 232: 69-76. 

Nelson, K.E., Weinel, C, Paulsen, IT., Dodson, R.J., Hilbert, 
H., Martins dos Santos, V.A., etal. (2002) Complete 
genome sequence and comparative analysis of the meta- 
bolically versatile Pseudomonas putida KT2440. Environ 
Microbiol 4: 799-808. 



Neumann, G., Kabelitz, N., Zehnsdorf, A., Miltner, A., 
Lippold, H., Meyer, D., et al. (2005) Prediction of the adapt- 
ability of Pseudomonas putida DOT-T1E to a second 
phase of a solvent for economically sound two-phase 
biotransformations. Appl Environ Microbioll^ : 6606-661 2. 

Nijkamp, K., van Luijk, N., de Bont, J.A., and Wery, J. (2005) 
The solvent-tolerant Pseudomonas putida S12 as host for 
the production of cinnamic acid from glucose. Appl Micro- 
biol Biotechnol 69: 1 70-1 77. 

Octave, S., and Thomas, D. (2009) Biorefinery: toward an 
industrial metabolism. Biochimie 91: 659-664. 

Palleroni, N.J. (2010) The Pseudomonas story. Environ 
Microbiol 12: 1377-1383. 

Pini, C.V., Bernal, P., Godoy, P., Ramos, J.L., and Segura, A. 
(2009) Cyclopropane fatty acids are involved in organic 
solvent tolerance but not in acid stress resistance in 
Pseudomonas putida DOT-T1 E. Microb Biotechnol '2: 253- 
261. 

Pinkart, H.C., Wolfram, J.W., Rogers, R., and White, D.C. 
(1996) Cell envelope changes in solvent-tolerant and 
solvent-sensitive Pseudomonas putida strains following 
exposure to oxylene. Appl Environ Microbiol 62: 1129- 
1132. 

Ramos, J.L., Duque, E., Huertas, M.J., and Haidour, A. 
(1995) Isolation and expansion of the catabolic potential of 
a Pseudomonas putida strain able to grow in the presence 
of high concentrations of aromatic hydrocarbons. J Bacte- 
riol 177: 3911-3916. 

Ramos, J.L, Duque, E., Godoy, P., and Segura, A. (1998) 
Efflux pumps involved in toluene tolerance in Pseu- 
domonas putida DOT-T1 E. J Bacteriol 180: 3323-3329. 

Ramos, J.L., Duque, E., Gallegos, M.T., Godoy, P., Ramos- 
Gonzalez, M.I., Rojas, A., etal. (2002) Mechanisms of 
solvent tolerance in gram-negative bacteria. Annu Rev 
Microbiol 56: 743-768. 

Ramos, J.-L., Marques, S., van Dillewijn, P., Espinosa-Urgel, 
M., Segura, A., Duque, E., etal. (2011) Laboratory 
research aimed at closing the gaps in microbial bioreme- 
diation. Trends Biotechnol 29: 642-647. 

Ramos-Gonzalez, M.I., Godoy, P., Alaminos, M., Ben-Bassat, 
A., and Ramos, J.L. (2001) Physiological characteri- 
zation of Pseudomonas putida DOT-T1E tolerance to 
p-hydroxybenzoate. Appl Environ Microbiol 67: 4338- 
4341. 

Revelles, O., Wittich, R.-M., and Ramos, J.L. (2007) Identifi- 
cation of the initial steps in D-lysine catabolism in Pseu- 
domonas putida. J Bacteriol 189: 2787-2792. 

Rodriguez-Herva, J.J., Garcia, V., Hurtado, A., Segura, A., 
and Ramos, J.L. (2007) The ttgGHI solvent efflux pump 
operon of Pseudomonas putida DOT-T1 E is located on a 
large self-transmissible plasmid. Environ Microbiol 9: 
1550-1561. 

Rojas, A., Duque, E., Mosqueda, G., Golden, G., Hurtado, A., 
Ramos, J.L., and Segura, A. (2001 ) Three efflux pumps are 
required to provide efficient tolerance to toluene in 
Pseudomonas putida DOT-T1 E. J Bacteriol 1 83: 3967- 
3973. 

Rojas, A., Duque, E., Schmid, A., Hurtado, A., Ramos, J.L., 
and Segura, A. (2004) Biotransformation in double-phase 
systems: physiological responses of Pseudomonas putida 
DOT-T1 E to a double phase made of aliphatic alcohols and 



© 2013 The Authors. Microbial Biotechnology published by John Wiley & Sons Ltd and Society for Applied Microbiology, Microbial 
Biotechnology, 6, 598-611 



610 Z Udaondo et al. 



biosynthesis of substituted catechols. Appl Environ Micro- 
biol 70: 3637-3643. 

Salzberg, S.L, Delcher, A.L, Kasif, S., and White, O. (1998) 
Microbial gene identification using interpolated Markov 
models. Nucleic Acids Res 26: 544-548. 

Sashidhar, B., and Podile, A.R. (2009) Transgenic expres- 
sion of glucose dehydrogenase in Azotobacter vinelandii 
enhances mineral phosphate solubilization and growth of 
sorghum seedlings. Microb Biotechnol 2: 521-529. 

Segura, A., Rojas, A., Hurtado, A., Huertas, M.J., and 
Ramos, J.L. (2003) Comparative genomic analysis of 
solvent extrusion pumps in Pseudomonas strains exhibit- 
ing different degrees of solvent tolerance. Extremophiles 7: 
371-376. 

Segura, A., Godoy, R, van Dillewijn, P., Hurtado, A., Arroyo, 
N., Santacruz, S., and Ramos, J.L. (2005) Proteomic 
analysis reveals the participation of energy- and stress 
related proteins in the response of Pseudomonas putida 
DOT-T1 E to toluene. J Bacteriol 187: 5937-5945. 

Segura, A., Rodriguez-Conde, S., Ramos, C, and Ramos, 
J.L. (2009a) Bacterial responses and interactions with 
plants during rhizoremediation. Microb Biotechnol 2: 452- 
464. 

Segura, A., de Wit, P., and Preston, G.M. (2009b) Life of 
microbes that interact with plants. Microb Biotechnol 2: 
412-415. 

Sikkema, J., de Bont, J.A.M., and Poolman, B. (1995) Mecha- 
nisms of membrane toxicity of hydrocarbons. Microbiol Rev 
59: 201-222. 

Surwase, S.N., and Jadhav, J. P. (2011) Bioconversion of 
L-tyrosine to L-DOPA by a novel bacterium Bacillus sp. 
JPJ. Amino Acids 41: 495-506. 

Tao, R, Tang, H., Gai, Z., Su, R, Wang, X., He, X., and Xu, P. 
(201 1 ) Genome sequence of Pseudomonas putida Idaho, a 
unique organic-solvent-tolerant bacterium. J Bacteriol 193: 
7011-7012. 

Tatusov, R.L., Fedorova, N.D., Jackson, J.D., Jacobs, A.R., 
Kiryutin, B., Koonin, E.V., et al. (2003) The COG database: 
an updated version includes eukaryotes. BMC Bioinformat- 
ics 4: 41 . 

Udaondo, Z., Duque, E., Fernandez, M., Molina, L., de la 
Torre, J., Bernal, P., etal. (2012) Analysis of solvent toler- 
ance in Pseudomonas putida DOT-T1E based on its 
genome sequence and a collection of mutants. FEBS Lett 
586: 2932-2938. 

Ussery, D.W., Kiil, K., Lagesen, K., Sicheritz-Ponten, T, 
Bohlin, J., and Wassenaar, T.M. (2009) The genus Burkhol- 
deria: analysis of 56 genomic sequences. Genome Dyn 6: 
140-157. 

Van Dillewijn, P., Caballero, A., Paz, J. A., Gonzalez-Perez, 
M.M., Oliva, J.M., and Ramos, J.L. (2007) Bioremediation 
of 2,4,6-trinitrotoluene under field conditions. Environ Sci 
TechnolM: 1378-1383. 

Volkers, R.J.M., de Jong, A.L, Hulst, A.G., van Baar, B.L.M., 
de Bont, J.A.M., and Wery, J. (2006) Chemostat-based 
proteomic analysis of toluene-affected Pseudomonas 
putida S12. Environ Microbiol 8: 1674-1679. 

Weber, F.J., Ooijkaas, L.P., Schemen, R.M., Hartmans, S., 
and de Bont, J.A. (1993) Adaptation of Pseudomonas 
putida S12 to high concentration of styrene and other 
organic solvents. Appl Environ Microbiol 59: 3502-3504. 



Wierckx, N.J., Ballerstedt, H., de Bont, J.A., and Wery, J. 
(2005) Engineering of solvent-tolerant Pseudomonas 
putida S12 for bioproduction of phenol from glucose. Appl 
Environ Microbiol 71: 8221-8227. 

Wierckx, N.J., Ballerstedt, H., de Bont, J.A., de Winde, J.H., 
Ruijssenaars, H.J., and Wery, J. (2008) Transcriptome 
analysis of a phenol-producing Pseudomonas putida S12 
construct: genetic and physiological basis for improved 
production. J Bacteriol 190: 2822-2830. 

Wu, X., Monchy, S., Taghavi, S., Zhu, W., Ramos, J., and van 
der Lelie, D. (2011) Comparative genomics and functional 
analysis of niche-specific adaptation in Pseudomonas 
putida. FEMS Microbiol Rev 35: 299-323. 



Supporting information 

Additional Supporting Information may be found in the 
online version of this article at the publisher's web-site: 

TextSI. Details on enzyme screens. 
Fig. S1. Catabolism of glucose as an example of sugar 
metabolism based on annotated genes. Glucose can be oxi- 
dized in the periplasm to yield gluconate or ketogluconate or 
in the cytoplasm upon phosphorylation. All three pathways 
converge at the level of 6-phosphogluconate that is metabo- 
lized via the Entner-Doudoroff pathway that yields chemicals 
that feed the Krebs cycle. The genes are annotated according 
to del Castillo etal. (2007). 

Fig. S2. Catabolism of L- and D-lysine as an example of an 
amino acid being used as a source of carbon through multiple 
confluent pathways. The set of genes and most of the inter- 
mediates were analysed in detail by Revelles and colleagues 
(2007) for the KT2440 strain. All genes shown in this pathway 
for DOT-T1 E are based on the identification of the homolo- 
gous genes. 

Fig. S3. Catabolism of lactic acid by P. putida DOT-T1 E as 
an example of the metabolism of an organic acid. Information 
was derived from Nelson and colleagues (2002). Lactate 
enables the growth of T1E in minimal medium with doubling 
times in the range of 120 ± 10 min. Inactivation of the HdD 
gene blocks the use of lactate as a C source. 
Fig. S4. Catabolism of aromatic compounds to a set of 
central intermediates. For the described reaction the corre- 
sponding enzyme(s) were identified as homologous to those 
described in detail by Jimenez and colleagues (2002) for 
KT2440 strains. A set of peripheral enzymes lead to the 
formation of a number of catechol-related compounds, 
which upon ortho or mefa-cleavage yielded Krebs cycle 
intermediates. 

Fig. S5. Genetic organization of the genes for the protocate- 
chuete central catabolic pathway. The genes were identified 
by blast analysis and the operon structure of the genes 
deduced from the overlapping nature of all the genes, except 
pcaF and pcaT; however, these genes form an operon based 
on RT-PCR analysis with RNA isolated from strains grown in 
4-hydroxybenzoate (unpublished). 

Fig. S6. Genetic organization of the duplicated cat genes 
located in two different chromosomal regions. One of the car 
clusters includes the catR gene that probably controls the 
expression of the two catABC operons based on high 
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sequence conservation of the promoter region in the two catB 
genes. 

Fig. S7. Genes and enzymes for the utilization of methionine 
as an S source by DOT-T1E. The gene products involved in 
the catabolism of methionine were identified based on blast 
analyses. 

Fig. S8. Set of biotransformation reactions to achieve the 
synthesis of adipate/3-oxoadipyl-CoA based on gene content 
in DOT-T1E strain. Multiple steps are represented by X 1 , X 2 
and X 3 . The loci are: 

X 1 . - T1 E-031 4, T1 E-0333, T1 E-0690, T1 E-0898, T1 E-1 004, 
T1E-1466, T1E1534, T1E-2226, T1E-2233, T1E-2397, 



T1E-3843, T1E-3853, T1E-3860, T1E-4416, T1E- 

5253, T1E-5478; X 2 . - T1E-0103, T1 E-1 864, T1 E-1 899, 

T1 E-2485, T1 E-2953, T1 E-3631 , T1 E-471 6, T1 E-4861 , T1 E- 

4862; X 3 . - T1 E-0369, T1 E-1 262, T1E-4182, T1E-4445. 

Table S1 . Pseudomonas putida DOT-T1 E enzymes of 

according to four digits EC classification. 

Table S2. C, N and S sources identified in silico and used by 

P. putida DOT-T1 E confirmed experimentally. 

Table S3. Dioxygenases encoded by the P. putida DOT-T1 E 

genome. 

Table S4. Monooxygenases encoded by the P. putida DOT- 
T1 E genome. 
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